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Abstract 

Two relativistic models for collapsing spheres at different stages of 
evolution, which include pre-relaxation processes, are presented. The 
influence of relaxation time on the outcome of evolution in both cases 

O . 

q-i is exhibited and established. It is shown that relaxation processes 



can drastically change the final state of the collapsing system. In 
particular, there are cases in which the value of the relaxation time 



^ j determines the bounce or the collapse of the sphere. 
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1 Introduction 



In the study of gravitational collapse, where dissipative processes play a fun- 
damental role, thermal relaxation time (r) is usually neglected. The reason 
for this may be found in the well known fact that, for most materials at 
usual laboratory conditions, r is very small as compared with typical scales 
of time of most self-gravitating systems. Thus, it is of the order of 10 -11 s for 
phonon-electron interaction and of the order of 10~ 13 s for phonon-phonon 
and free electron interaction at room temperature Q. 

There are however situations where relaxation time may not be negligible. 
Thus, for example, for superfluid Helium, r is of the order of 10 -3 s for a 
temperature 1.2 K ||. 

Also, and much more important in the context of the problem considered 
here, in cores of evolved stars the electron gas is highly degenerate. Then, 
since the quantum cells of phase space are filled up, such that collisions in 
which the momentum is changed become rather improbable, the mean free 
path increase considerably, increasing r thereby. 

Thus, for example, for a completely degenerate core of radius ~ 10 _2 i? Q 
at T ~ 10 7 K, the thermal relaxation time is of the order of one second [[]]. 
However, the order of magnitude of relaxation time used here is much lower 
than that. 

Using the expression for the thermal conductivity by electrons in neutron 
star matter introduced by Flowers and Itho |BL we found (see below) 
without difficulties relaxation times in the range [10 -6 s, 10~ 4 s]. 

There have been recent calculations on collapsing systems which include 
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pre-relaxation processes [7]]. However, either they consider a constant 
value for the conductivity || or they are calculated for a single value of r [[?]], 
excluding the possibility of comparing the evolution for different relaxation 
times, and thereby, of assessing the influence of r on the evolution of the 
collapsing object. 

In this paper we describe the evolution of two different self-gravitating 
systems dissipating energy through a radial heat flow vector. 

Modeling is achieved by using the HJR formalism || |9], lOfl . However, 



unlike Di Prisco, Herrera and Esculpi || we shall not consider a constant 
value for thermal conductivity and the astrophysical setting is much more 
realistic. Also, unlike Martinez @ we shall follow the approach introduced 
by Di Prisco et al || which allows for comparing the evolution for different 
values of r. 

The first model is based on the well known Tolman VI solution and 



may be accomodated to describe to some extent, the core implosion and 
subsequent bounce, at earlier stages of a supernova explosion. 

The second one is more suitable for describing the Kelvin- Helmholtz phase 
of the birth of a neutron star [Q. However, it is important to emphasize 
that in spite of the fact that the order of magnitude of different physical 
variables are well within the range of expected values, our main goal here is 
not to present a detailed modeling of those scenarios, but to bring out the 
relevance of pre-relaxation processes in situations when degenerate cores are 
involved. 

The paper is organized as follows. In the next section the field equations 
and conventions are presented and also a brief resume of the HJR formalism 



is given. The heat conduction equation is presented in section 3, as well as 
the approach to operate the HJR formalism. 

Section 4 is devoted to the description of the two models considered in 
this paper. 

Finally, some conclusions and comments are included in the last section. 

Partial differentiation with respect to u and r are indicated by means of 
the subscripts and 1 respectively. The subscript a denote that the quantity 
is evaluated at the surface. 



2 The HJR formalism 

We shall consider a non-static distribution of matter which is spherically 
symmetric, and consists of fluid (which may be locally anisotropic) and heat 
flow. In Bondi coordinates O, ITU , the metric takes the form 



ds 2 = e 2ii 



V 2 

— du + Idudr 
r 



r 2 



(dO 2 + sin 2 6d<f) 2 ) , (1) 



where u = x° is a time like coordinate (g uu > 0) , r = x 1 is the 
null coordinate (g rr = 0) and 9 = x 2 and <fi = x 3 are the usual angle 
coordinates. The ^-coordinate corresponds to the retarded time in flat space- 
time and the metric functions (3 and V are functions of u and r. A function 
m(u,r) can be defined by 

V = e 2f3 (r -2m(u,r)), (2) 
which is the generalization, inside the distribution, of the "mass aspect" de- 



fined by Bondi et al [13]. In the static limit it coincides with the Schwarzschild 



mass. On the other hand, the radiation coordinates u, r, 9 and 4> can be re- 
lated to Schwarzschild ones (T, R, 6, <£>) by 

T = u+ [ r ?-dr, (3) 
Jo V 

R = r, = 6, $ = 0, (4) 
and to local Minkowskian coordinates (t, x, y, z) by 

dt = e^^-du + sj^dr), (5) 

dx = e^J^dr, (6) 
dy = rd6, (7) 
dz = rsin9d(p. (8) 
At the outside of the fluid distribution the space-time is described by 



the Vaidya metric ||14| , which in Bondi coordinates is given by (3 = and 
V = r — 2m(u). 

For the matter distribution considered here, the stress-energy tensor can 
be written as j7| 

= (p + P±)U ll U u - P L g^ v + (P r - P±)x»Xu + 2Q^U v) , (9) 

where p, P r , P± are the energy density, radial pressure and tangential pres- 
sure respectively as measured by a Minkowskian observer in the Lagrangean 
frame. The stress-energy tensor outlined above is obtained as measured by 
an observer using Bondi coordinates. For this observer the four- velocity U^, 
and the heat flux vector are given by 



and 



Q 1 = Qe 



-0 



V VI 



CO 



id 



2 iv 1 

7(1 _ ^2)1/2 > ' 



whereas — ~Q^/ Q, Q = \J Q^Q^ is the heat flow, and oj is the velocity 
of matter as measured by the Minkowski observer defined by 

The Einstein field equations, inside the matter distribution, can be writ- 



ten as |T(J: 



47rr(r — 2m) 



(-m e~ 2 P + (1 - 2m/r)m^\ = —^(p + 2uQ + P r u 2 ), (12) 

^=ih (p - Q{l - u) - p - u) - (13) 
^=^- 2G+p ^ (14) 

_ ?Bfl + ' (1 - 2*) (2ft, + 4# - ^ ) + 3ft " " f " ft " = Fx, (15) 
47T 57T r r 87rr 

while outside matter, the stress-energy tensor takes the form 



?V = -^m„^: (16) 



and the only non trivial Einstein equation reads 

2 / 2m(w)\ 
m = -47rr 2 £: I 1 — J . (17) 

where e is proportional to the energy density of the radiation traveling in 
the radial direction - see Bondi |L5| for details, and subscripts 0,1 denote 
derivatives with respect to u and r respectively. 
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In the HJR formalism R 0, [TO] , one introduces the concept of effective 



energy density and effective pressure. From (|i~3"D, the mass function can be 
expressed as 

rr 

rh= Anr 2 pdr, (18) 
Jo 

where 

p=-J—(p-Q(l- u )-P rU ) i (19) 
l + u 

is the effective energy density, which in the static limit reduces to the energy 
density of the system. 



or 



In a similar way, we can define the effective pressure. From (|14]) one has 

r r 97rr 2 1 1,1 



» = Ih^ + * (21) 



with 



P = -—{-u P -Q{l-u)+P r ) (22) 

1 + UJ 



being the effective pressure. This one, as the effective energy density, only 
has a clear physical meaning in the static case, in which it reduces to the 
radial pressure. 

Matching the Vaidya metric to the Bondi metric at the surface (r = a) 
of the fluid distribution implies f3 a = f3(u, r = a) = and the continuity of 
the mass function fh(u,r). The continuity of the second fundamental form 
must be demanded as well, leading to condition 

1 _ 2 ^« > \^3__ (23) 

O / Pa + Pa 



- see Herrera and Jimenez jlB] for details - where overdot denotes derivative 
with respect to u. 

The well-known condition 

Qa = Pra, (24) 

for radiative spheres |l7j can be obtained from the coordinate transformation 



(||). Effectively, the velocity of matter in Bondi coordinates can be written 

as 

t = -r^' < 25 > 

du r 1 — uj 

evaluating the last expression at the surface and comparing it with fl23|) it 
follows 

Pa = -UJaPa- (26) 

It is easy to show, using ( |TT?D and (^), that this condition is equivalent to 
(II- 

The HJR method is based in a system of three surface equations which 
allows us to find the u dependence of the functions f3 and in present in the 
field equations ( p^ - |T5|) . 



To derive the surface equations, usually five dimensionless variables are 
defined 

a ™ u 2M n 1 

A = t^t^ M = —- u = — - F = 1 VI = , 27 

m(0) m(0) m(0) A 1 - uo a y ' 



where m(0) is the initial total mass of the system. Then, using fl2?| ) and (|26|). 
( |2T)| ) yields the first surface equation. Using the functions just defined into 
(|23|) we get the first surface equation 

A = F(n-l), (28) 



which gives the evolution of the radius of the star. From now on, a dot 
denotes derivative with respect to the dimensionless variable u/m(0). 

The second surface equation can be obtained from the luminosity evalu- 
ated at the surface of the system. The luminosity as measured by an observer 
at rest at infinity reads 

L — —M — - E X9 =EF = E(2Q - 1)F = 4nA 2 Q a (2fi - 1) F, (29) 

(i + z a y 

where z a refers to the boundary gravitational redshift, E is the luminosity 
as seen by a comoving observer, and E is the luminosity measured by a non 
comoving observer located on the surface. Using relationship (|T2]) and (^9]) 
together with the first surface equation we obtain the second one as 

2L + F(l-F)(n-1) , , 

F = L, (30) 

which expresses the evolution of the redshift at the surface. 

The third surface equation is model dependent. For anisotropic fluids the 
relationship (T^) a = can be written as 

— + — ^ + Fn 2 ^- — Ftt— = 

F il p a p a A p a 

il 2 A p a Ap a 



being 



Expression (|3l|) generalizes the Tolman-Oppenheimer-Volkov equation to the 
non-static radiative anisotropic case. 
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The HJR method M allows us to find non static solutions of the Einstein 
equations from static ones. The algorithm, extended for anisotropic fluids, 
can be summarized as follows 0, 3 

1. Take a static but otherwise arbitrary interior solution of the Einstein 
equations for a spherically symmetric fluid distribution 



2. The effective quantities p = p(u,r) and P = P(u,r) must coincide 
with p st and P st respectively in the static limit. We assume that the in- 
dependence in effective quantities is the same that in its corresponding 
static ones. Nevertheless, note that junction conditions in terms of 
effective variables, read as (0). This condition allows us to find out 
the relation between the w-dependence of p = p(u, r) and P = P(u, r). 

The rationale behind the assumption on the r dependence of the ef- 
fective variables P and p, can be grasped in terms of the characteristic 
times for different processes involved in a collapse scenario. If the hy- 
drostatic time scale Thydr, which is of the order ~ l/y/Up (where G 
is the gravitational constant and p denotes the mean density) is much 
smaller than the Kelvin- Helmholtz time scale ( Tkh ) , then in a first ap- 
proximation the inertial terms in the equation of motion can be ignored. 
Therefore in this first approximation (quasi- stationary approximation) 
the r dependence of P and p are the same as in the static solution. 
Then the assumption that the effective variables have the same r de- 
pendence as the physical variables of the static situation, represents a 




Pst 




(33) 
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correction to that approximation. 

3. Introduce p(u,r) and P(u,r), into (|18|) and ([2l|) to determine m and 
(3 up to three unknown functions of time. 

4. The three surface equations form a system of first order ordinary dif- 
ferential equations, by solving it we find the evolution of the radius, 
A(u), and two unknown functions of time. These ones can be related 
with the u-dependence of p = p(u, r) and P = P(u, r). 

5. There are four unknown functions of time (A, F, Q and L). Thus, it is 
necessary to impose the evolution of one of them to solve the system 
of three surface equations. Usually the luminosity is taken as an input 
function because it can be found from observational data. 

6. Once these three functions are known, it is easy to find m and j3. 
Therefore, the interior metric is completely defined. 

7. Now, the left-hand side of the Einstein equations (|l^)-(|l5|) is known. 
However, the right-hand side of these equations contain five unknown 
quantities (u, p, P r , P± and Q). Thus, it is necessary to supply another 
equation to close the system of field equations. In the anisotropic static 



case the equation of hydrostatic equilibrium reads [18 



2 V 2 J \ r-2m 

This expression is usually generalized, in the context of HJR method, to 
non-static cases by substituting the physical quantities by the effective 
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variables H, IT9 



Pi 



' p + P\ (m + Anr^P^ 



r — 2m 



(35) 



Now, the Einstein equations, supplemented with (p5[), form a closed 
system of equations, and quantities u, p, P r , P± and Q can be found. 



3 System of Surface Equations: Heat con- 
duction equation. 

As we mentioned in the previous section, equations (H), (|30|) , and con- 
stitute a system of surface equations. However, if we desire to study the 
influence of the relaxation processes during the collapse, it is necessary to in- 
troduce an hyperbolic transport equation. Recently, Di Prisco, Herrera and 
Esculpi || have shown the importance of the relaxation processes involving 
the heat flow using a Schwarzschild based model with constant thermal con- 
ductivity. Our aim is to apply the same procedure to different and more 
realistic models in order to discern the effects that are model independent 
and how sensitive is the final state to the presence of relaxation processes. 

The system of surface equations can be solved for a given luminosity pro- 
file. Thus, the temperature could be found if r were known 0. Nevertheless, 
this procedure does not give us information about the influence of the relax- 
ation processes on the luminosity profile. Because of this, it seems more 
convenient to our purposes to follow the method outlined by Di Prisco et al 

I- 
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First of all, we assume the evolution of the heat flow governed by the 
Maxwell-Cattaneo transport equation, 

JU V ' 



dQ v 

as 



T„ - T- 



ds 



(36) 



where 



dQ v 
ds 

dU» 
ds 



and x is t ne thermal conductivity coefficient. The non vanishing covariant 
components of the four acceleration read: 



dU r 
ds 



1 



1+LJ 



--23 - l ~ 2 ™ 1 
2r 2(r - 2m) 



+ re 



-23 



1-U 



m 



1 



uuJi + re 



-23 



1 + uj) (r — 2m) 2 
u 



r — 2m 



and 



dU u 



2/3 / _ 2m\ 



2m \ uj dU r 



(37) 



(38) 



ds \ r J 1 — to ds 

It is easy to demonstrate that in our case, eq.(|36|) has only one independent 
component. Thus, evaluating one of them at the surface we obtain 



r Q a +Q a JF(2n - 1 



Xa 



T a -T la F(2n - 1) - T a F{2Vt - 1) + F 2 A(2Vt — 1) + 

(39) 



where we have used (gTj), (0), (0), and (0). 
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We shall assume that the thermal conduction is dominated by electrons. 
Thus, the thermal conductivity is given by expression |4], U 



X ~ 10 23 — erg. s" 1 cm" 1 K" 

^8 



(40) 



where the energy density and the temperature are given in 10 14 g cm -3 and 
10 8 K units respectively. In the HJR formalism the initial mass is normalized 
to unity (^7]). Therefore, all quantities that, in geometrized units, have 
dimensions of a power of length, are within this framework dimensionless. 
The expression ( |4"0| ) for the thermal conductivity coefficient can be rewritten 
as 

X*%L (41) 



being C = 2.5126 x 10~ 20 , p the dimensionless energy density, and T the 
temperature in Kelvin. The initial mass has been written in terms of the 
solar mass as M Q = £M Q , where £ is a numerical factor. 

The energy density in the surface can be written in terms of the heat flow 
and the effective energy by means of the definition ([191) . Evaluating it at the 
surface and using ( p7|) and junction condition (|24| ) we obtain, 

■2ft - r 



Pa= Qa + Pa 



ft 



(42) 



Then, evaluating (JI) at the surface, the heat transport equation at the 



surface (39) reads 



r Q a +Q a JF(2n-l 



^-^F(2ft-1)-F(2ft-1) 



C 



l-F 



Qa + Pa 
ft 



2ft - 1 



ft 

+ 



2AF J 2ft - 1 F 2 A{2VL - 1 



+ 



ft 



F(2ft 
(43) 
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3.1 Boundary condition for the transport equation 



Expression (|43|) relates the heat flow and the temperature at the surface to 
the temperature gradient through quantities present in the surface equations 



(HD, ©, and 

The connection, at the surface, between the heat flow and the temperature 
can be found by means of the effective temperature T e ff. This one is defined 
by means of the luminosity perceived by an observer located momentarily on 
the surface, 



E 



47TTVZJ, 



(44) 



being a the Steffan-Boltzmann radiation constant. Thus, the effective tem- 
perature would be the temperature at the surface of the star if it would 
radiate as a black body. This concept is linked up with the theory of stellar 



atmospheres by means of the photosphere [PQj p. 70], p. 586] and |22 



p. 295]. The photosphere comprises the most external layers of the star, and 
its thickness is determined by the optical depth for photons. At the surface 
this optical depth vanishes and the effective temperature is related to the 
material one by the expression 

1, 



rpA 

a 



2 ef f 



Then, by substitution of ( f4"5|) and ( pH] ) into ( p9|) we obtain 

2a 



Qa 



rpA 

a 2tl 



1 



(45) 



(46) 



Introducing this expression in (|43|) , we obtain the fourth surface equation 



4 £-2 



2Q 



At 



Pa 



C /2Q - 



2„.5 



4r VI 



(2tt 



2\JF 
v / 2TT= 



15 



c 
1 

where 



2y 4 £ 2 „ (20. -I 



20 - 1 + Pa V o 



(2/^(20-1)+^), (47) 



*=(^rJn + W« 2 + (48) 



the function y has been defined as 

y 4 = Ctf, (49) 

and the constant ( = crM| ~ 3.4097 x 10~ 54 K" 4 is calculated taken M 
in geometrized units. The luminosity L, present in (f|3), can be written in 
terms of y by means of fl2~9|) , d45f) and (|4]) 

L = 8ttA 2 ?/ 4 ^ 2 F. (50) 

Note that all terms in ( fl"7l) are dimensionless, including the relaxation time 
t, which is related to the physical r(= r^) by means of 

T = T ph -^- ~ 2.0298 x 10 5 ^ s~\ (51) 

where c and G denotes the speed of light and the gravitational constant, and 
m is the total mass in grams. 



3.2 The DHE approach. 

As we mentioned above we shall use the approach introduced by Di Prisco et 
al (DHE approach henceforth) to solve the system of surface equations. 
This system of equations usually is composed of three differential equations 
which give the evolution of three quantities: The radius (|28"D, the boundary 
gravitational redshift ( [30]) and the collapse velocity of the surface (|3T|). In 
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order to study the effect of the thermal relaxation processes on the system it 
is necessary to include an additional equation accounting for the evolution of 



the temperature at the surface. Thus, expression ([47]) constitutes the fourth 



surface equation. The method can be described as follows: 

1. Solve, as usual, the system of surface equations (|28|) , (|30|) , and (|3lD for 
a given luminosity L. 



2. Introduce the results into the fourth surface equation (|47|), and solve 
if for t = (Eckart's law). Since L is known, equation fl50|) allows us 
to find the temperature evolution and equation ( JT7D provides us the 
temperature gradient. 

3. Use the evolution found for T\ a in the resolution of the complete sys- 
tem of surface equations (|28l), (|3~0D , (|3~T1), and ( |47|) for non vanishing r 



(Maxwell-Cattaneo law). So that the evolution of A, F and L is found 
for a given value of r and we are able to compare it with that found 
for t = 0. 

The method described above must be self consistent, i.e. we must recover 
the luminosity profile assumed in the first point if a vanishing relaxation time 
is imposed in the third one. 

The procedure adopted here assumes that the temperature gradient is 
independent of the adopted value for r. This hypothesis can be justified 
in the following terms: The existence of a temperature gradient is caused 
by local processes. As a consequence of its presence a non-vanishing flux 
appears, giving rise to dissipative processes. The relation between cause and 
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effect is, in this case, non linear due to the fact that dissipative processes 
can change the local situation. On the other hand, relaxation processes 
act directly on dissipative processes. Thus, they must affect in a indirectly 
way to the local situation. Nevertheless, the effect on local processes must 
be less important than on dissipative ones having, then, a little influence on 
temperature gradients. In first approximation we assume that, for time-scales 
of the order of r, the influence of the relaxation processes on temperature 
gradients can be neglected. 

In order to ascertain the goodness of this assumption, we have numerically 
integrated the system of surface equations Q28D, Q50|), ([H]) and (0), for the 
models presented below, assuming the same luminosity profile for a wide 
range of values of r. The differences found among temperature gradients for 
different r can be neglected. This fact reinforces the adopted hypothesis. 
The unsensibility of the inner temperature to the relaxation time for a fixed 
luminosity [|7| points in the same direction as well. 



4 The models 



The HJR method starts from a static solution of the Einstein equations (sec. 
II). We shall study the evolution of two different models. The first one, 



corresponds to the well-known Tolman VI solution |]TTJ , while the second one 



comes from a static solution due to Gokhroo and Mehra 23 
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4.1 Tolman VI- type solution 



The equation of state obtained from the Tolman VI static solution ap- 
proaches, in the core of the star, to a highly compressed Fermi gas. The 
energy density and radial pressure are 

3 



p(r] 



St 



and 



P(r 



P 



567rr 5 



9a - 9r 



St 



9a — r 



(52) 



(53) 



respectively. In order to generate non static solutions, Tolman VI model has 
been used often in the HJR framework as a departing solution ||. Following 
the HJR method, the effective quantities are taken as 

3G(u) 



and 



P 



P= P - 
3 



9K(u)r 



1 - K{u)r 

Using expressions fl54|), ( |I8D and fl2~7|) the function G{u) reads, 

1 -F 



G 



24tt 



(54) 



(55) 



(56) 



Function K(u) follows, as a function of Q, from junction condition (p6|) , (|5^) , 
(H) and (|3), 

By means of the last four expressions, the system of surface equations 
(M), (M), (P), and (O) reads 



A = F(n- 1), 
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(58) 



2L + F(1- F)(fl-1) 



and 



A 

ft _ F(l - F)(4fi - 3) (4ft - 1) 
' ft ~ 8A 



Am 



;i - FY 

2AVL (2ft -i)A' 



(59) 
(60) 



4t2 



2ft- 1 



4r 



C 
1 



I) 



+ 



C n-F\ /2ft- 1 
f V 8ttA 2 ) V ft 

1 - F\ /2ft - 1 



4r ft 



2VF 



(2ft - 1) 2 v^TT^T 
(2/^(20 -l)+y$), (61) 



2ft - 1 ' V 8vr/L 2 y v ft 
This model can be applied in the context of the bouncing of the core 
during a supernova explosion. It is important to note that the value of the 
quantity M/A plays a fundamental role in the evolution of the system. In 
particular, for the model under consideration the bounce of the surface occurs 
if M/A < 3/14 ||. As Di Prisco et al || showed, this quantity is influenced 
by the inclusion of relaxation processes. If this effect is model independent 
then, in some cases, it will cause a drastic change in the final state of the 
collapsing core. 

In order to evaluate the influence of non vanishing r in the evolution of 
the system, we assume as initial configuration: 



A = 6, ft D = 0.878497, F = -, and M D = 1.3M . 

3 



(62) 



This corresponds to a core with an energy density about 10 14 g cm -3 in the 
external layers, whereas its radius is 11521 m. 

In the case of vanishing relaxation time, the luminosity is taken as 



2M 



s/. 



r =Vte-^ 



X°7T 



(63) 
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where M r is the total radiated mass, and x is an arbitrary parameter which 
determines the decaying time of L. In this has been taken as 10 and 

M r = 8.407 x 10 -6 . This brings a maximum luminosity around 1.5 x 10 53 
erg s _1 , and a characteristic time for L close to 0.5 ms. The total radiated 
energy is of the order of 10 49 erg. 

Following the DHE method described in the previous section, the system 
of surface equations (|5B])-(l6T|) can be solved for r = using fl63|). This allows 
us to find the temperature gradient in the surface, which is used to solve 
dD-(|ID for r ^ 0. 

The luminosity profile is displayed in fig. 1 for several values of r. The 
maximum luminosity decreases if r grows. Thus, the higher r, the softer 
collapse. Other effect that can be inferred from figure 1 refers to the width 
of the pulse of luminosity. This one is larger for larger r's. The same influ- 
ence can be observed in the temperature of the surface (fig. 2). Note that 
the maximum luminosity and temperature takes place at larger times as r 
increases. This is due, as corresponds to a relaxation process, to the fact 
that for larger values of the relaxation time, the system requires more time 
to establish the heat flow. 

As we mentioned above, the quantity M/A is greatly influenced by the 
presence of relaxation processes (figure 3). An interesting effect must be 
noted in its evolution. The M/A profile presents a sort of bifurcation for 
T — T uf ~ 3.5 x 10~ 6 s. If r > Tuf, the core evolves towards a more compact 
final state than the initial one. The situation is reversed for r < r^f. Thus, 
small values of r have a great influence in the final state. Something similar 
occurs in the evolution of the radius (figure 4). The bounce of the surface 
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takes place if r < r^/, whereas for r > Tuf the collapse is ensured. This 
fact can be explained taking into account the influence of the luminosity in 
the evolution of the collapse. In the present model, and for r = 0, the pulse 
of luminosity is the responsible of the bounce of the surface. In case when 
luminosity is low, the explosion won't take place. The presence of relaxation 
processes imposes a lower luminosity. Therefore, for r > Tuf the luminosity 
will be insufficient to stop the collapse, and the surface won't bounce. There 
are cases in which, as in the studied model, small variations in the luminosity 
make the bounce impossible. 

The influence of thermal relaxation processes seems to be important in 
the evolution of the collapsing system. Nevertheless, the question is if the 
value found for ruj is close enough to the expected value for the relaxation 
time. Assuming that the thermal conduction is dominated by electrons, the 
relaxation time can be roughly estimated, in the limit of high frequencies, by 
the expression || 

10 20 

r ~ TV S ' (64) 
where T is measured in K and thermal signals travel with speed v given in cm 
s _1 . This expression is only valid assuming p and x as constants. However, 
it can be taken as an indicative of the order of magnitude for r. Assuming 
T ~ 10 10 K, and that thermal signals propagating approximately at 10 3 cm 
s _1 (the second sound speed in superfluid helium) we obtain r ~ 10~ 6 s. 
Thus, the obtained Tuf must be considered as a realistic value for r and the 
effects described above take a special importance. 

The total radiated mass, as a function of r, for two different initial and 
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boundary conditions is presented in figure 5. One of them (TOL06 hence- 
forth) corresponds to initial conditions ( |62"D and boundary condition (|63|). 
The other (TOL04 from now on) has been solved using 

A = 6.5, Q = 0.845215, F = — , and M = 1.3M , (65) 

13 

as initial conditions, and assuming a Gaussian pulse profile for the luminosity 
for vanishing r. The total radiated mass, if r = 0, is M r = 10~ 4 M o , and 
the maximum luminosity is about 7.2 x 10 54 erg s _1 . The point marked A 
in figure 5 indicates the limit of validity for TOL04 model. The Tolman 
VI solution fails for large times. In particular, TOL04 model fails before 
the luminosity may vanish if r is larger than ~ 0.07 ms. TOL04 model 
is interesting because it clearly shows the non linearity, between the total 
radiated mass, and r. This effect occurs in TOL06 model as well, but is only 
visible for very small relaxation times and the effect is not so accused. In 
spite of the maximum luminosity decreasing as r increases, the total time 
in which the system radiates grows with the relaxation time, and the total 
radiated mass too. Nevertheless, the maximum value of the total radiated 
mass is bounded. This upper limit to M r is reached if relaxation time is of 
the order of the width of the luminosity profile. 

4.2 Gokhroo &; Mehra-type solution 

This solution corresponds to an anisotropic fluid with variable energy density. 
The choice of this solution is based on the fact, that it leads to densities 
and pressures similar to the well-known Bethe-Borner-Sato (BBS) Newtonian 
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state equation |ZT , |24] , [25[ |. Its generalization to the HJR method has been 
done recently with good results 0. 



In this solution the energy density and radial pressure are assumed as [^3 

/ ,v,2\ 



p(r) st = p c \l-k— 



a' 



and 



2m\ ( r 
; (!-—) > „ 



(66) 



(67) 



where the constants k and n are within the range < k < 1 and n > 1. 
The central energy density p c , and central radial pressure P c are related by 
expression 

P c = Xp c . 



(68) 



The tangential pressure is given by 

r,„, { p + P r \ / m(r) + Airr 3 P r 
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The effective quantities are defined, following the second point of the HJR 
method, as - see Martinez |7j for details - 

/ ji \ 



P = Pc(u) 1 - K(u) — 



and 



a* 



P = P c (u) 1-2= )h-G(u) r - 2 



Where the functions p c (u) and P c (u), are defined by 



and 



Pc{u) = p, 



Priu) = P, 



K{u) 



K(u) 



K 



K(Q) '"K 
The functions of time K(u), and G(u) are found to be 
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The expression for tangential pressure can be found following the point seven 
of the HJR method. 

Thus, the system of surface equations for this model is 



A= F(fi-l) 



F=-[2L + F(l-F)(n-1)] 



(78) 
(79) 
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n= -In + ^~ 2i % + ^ 
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and $ is denned in (|48|). 

We have integrated the system of surface equations for this model with 
the initial data A(0) = 10; ft(0) = 1, and with n = 1; A = 1/3; if(0) = 0.999. 
Also, the inital mass is 1.3M . This corresponds to a star with an initial 
radius of 19.201 meters, a central density equal to 2.17 x 10 14 g cm -3 and 
a surface density equal to 2.17 x 10 ll gcm~ 3 . For the luminosity we have 
taken a gaussian pulse with center at u — 200 and width equal to 20 (in 
dimenssionless it-units), which is equivalent to a width of 0.13ms. 

In the case r = 0, the collapse spans over u ~ 400 (2.5ms). As r increases, 
that time also increases and may be as large as 32ms for r = 900 (5.7ms). 

In figure 6 we plot the ratio of the total radiated mass divided by the 
total radiated mass for r = 0, as function of r, for different radiated mass in 
the r = case. 



26 



The sensitivity of this ratio to different values of r is clearly exhibited. 
This fact is also present in the Tolman VI model discussed in the previous 
section (fig.5). However in this model non linear effects seem to be stronger 
than in the precedent one and the ratio may be larger or smaller than one, 
depending on the total radiated mass. 

The evolution of the radius is exhibited in figures 7 and 8 for different 
total radiated mass (1% and 0.001% of the total mass) and different r's. The 
dependence of the final value of the radius on r, is linear for low emission 
(0.001%) whereas it is not for the stronger case (1%). 

The ratio Mf ina i/Af ina i behaves essentially as the other surface variables 
(fig.9). The resulting object is more compact for larger luminosities, however 
as r increases, the final configuration is less compact and it may happen that 
for sufficiently large r's, the object may be at the end of it's evolution less 
compact than in the r = case. 

Finally, figures 10 and 11 show the evolution of luminosity profiles for 
different radiated mass and different r's. As in all known models, larger r's 
means longer emissions and more flattened pulses. This effect is sharper for 
stronger emissions. 

5 Conclusions 

It has been the purpose of this paper to exhibit the relevance of thermal 
relaxation time in the problem of collapse, by means of two models of ra- 
diating spheres. It is important to emphasize that this relaxation time is 
systematically neglected in collapse calculations, where processes may occur 
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on time scales which may be of the order of magnitude of (or at least not 
much larger than) relaxation time, leading thereby to incorrect conclusions. 

One of the models (Tolman VI) is suitable for describing the evolution of 
the core at the earliest stages of a supernova explosion. 

The second one (G-M) is more adapted to describe the Kelvin-Helmoltz 
phase in the birth of a neutron star. 

However, due to time restrictions in numerical calculation, in this model 
(G-M), we have integrated the surface equation over a period of about tens 
of milliseconds, instead of tens of seconds, which is the typical time of the 
Kelvin-Helmholtz phase in neutron star formation. Consequently all times 
have to be scaled by the same factor. 

Preliminary results show that in what the relaxation time concerns, its 
influence on the evolution of the object is not qualitatively changed by this 
"scaling" . 

In the Tolman VI model, the point to emphasize is the bifurcation in- 
troduced by changes in r of the order 1CT 6 s. Such values are within the 
range of possible values of r, and therefore the study of the evolution of such 
system seems to require a good account of pre-relaxation processes. 

In the second model the dependence of the final configuration on r may 
be non-linear and is affected by the total emission. 

At any rate the final configuration is clearly r-dependent. 

It is worth mentioning that the hydrostatic time scale for the second 
model is of the order of 0.29ms, which clearly indicates that in these cal- 
culations the hydrostatic approximation is not a very good one. This is an 
additional argument to use the HJR-formalism, which in some sense may be 
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envisaged as a "correction" to hydrostatic approximation (see point 2 of the 
algorithm in section 2). 

We would like to conclude with the following comment: In a collapse cal- 
culation, neutrino transport plays an important role, the reason to overlook 
this issue here, resides in the fact that we are not concerned with the problem 
of modeling gravitational collapse, but with the influence of thermal relax- 
ation time in the outcome of evolution. Conduction associated with trapped 
neutrinos is certainly of greatest relevance, and it should be very interesting 
to find out the role played by relaxation time of that process. This however 
is out of the scope of this paper. 
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Figure captions 

Figure 1.- Luminosity profile as a function of time. In all figures the values 
of the relaxation time are given in miliseconds. 

Figure 2.- Evolution of the temperature in the surface for different values 
of r. 

Figure 3.- Evolution of M/A for different values of r. Note the bifurcation 
between r = 0.0032 ms and r = 0.0038 ms. 

Figure 4.- Radius as a function of time. The bounce is only possible for 
t < 0.0035 ms. 

Figure 5.- Ratio between total radiated mass and total radiated mass for 
r = as a funtion of the relaxation time. Dashed and solid lines 
correspond to TOL04 and TOL06 models respectively. TOL04 model 
fails for r larger than the corresponding to point marked A. 

Figure 6.- Same as fig. 5 for G-M model. Curves are labeled with the 
percentage of radiated mass in the case r = 0. 

Figure 7.- Variation in the radius for M r = 10~ 5 (0.001%) 

Figure 8.- Same as fig. 7 for M r = 0.01 (1%). Note that in the low emission 
case (fig. 7) the variation is given in cm, while in this case AA is given 
in meters. 

Figure 9.- Ratio between final value of M/A and M/A for r = as a 
function of r. The labels in each curve means the same as in fig. 6. 
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Figure 10.- Luminosity profile for different values of r. The radiated mass 
in the case r = is a 0.001% of the initial mass. 

Figure 11.- Same as fig. 10. The radiated mass for r = is a 1% of the 
initial mass. 
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